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SUMMARY 


The conditions for similarity of flow for a compressible fluid with 
heat transfer are examined with particular reference to variation of 
stagnation temperature and scale of the model. Similar flows are pos- 
sible only when certain conditions are fulfilled by the intrinsic prop- 
erties of the gas . The restrictions on the characteristics of the free 
stream a.nd on the scale of the model necessary to obtain similar flows 
are pointed out. For similar flows in air, the following expression 
gives the ratio of local rates of heat transfer on the model to those on 
the prototype: 


^tunnel _ | ' °tunnel \ L f light 

flight ^flight/ L tunnel 

where 

H heat transfer per unit area per unit time 
T q absolute stagnation temperature 
L typical length 

Important simplification of the experimental procedure for tests simu- 
lating flight at high. Mach numbers is made possible by divorcing the 
fluid-heat -conduction problem from the heat-radiation problem and inter- 
preting the wind-tunnel results by means of similarity relations developed 
in this paper. 

INTRODUCTION 

The conditions for similarity of flow between model and prototype 
are of basic importance in problems associated with model testing. The 
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usual wind-tunnel test data are generally assumed -to "be strictly appli- 
cable to flight conditions when the values of the Reynolds number and 
Mach number are the same for the prototype in flight as for the model in 
the wind tunnel. Identity of the Reynolds number and Mach number ceases 
to be a sufficient condition for flow similarity if the effects of heat 
transfer are appreciable . 

The effects of heat transfer on the flow pattern can. be important 
at high Mach numbers. For example, calculations by Lees (reference l) 
have indicated that heat transfer from the boundary layer to the body 
has a marked stabilizing effect on_the l aminar boundary layer. The 
stagnation temperature corresponding to a Mach number of 4.0 in air and 
a free-stream temperature of 400° F absolute is l68o° F absolute, and the 
body surface will tend to assume a temperature close to this value. At 
such high temperatures, however, the amount of heat-dissipated per unit 
area from the body by radiation is large. Under steady conditions, a 
corresponding amount of heat must be transmitted to the - body from the 
surrounding fluid. Wind-tunnel measurements of the drag of bodies at 
high Mach numbers should, therefore, be made with rates of heat transfer 
to the surface corresponding to flight conditions in order for the tran- 
sition point to occur at the right place and, hence, for the skin friction 
to have the proper value . 

Exact duplication of flight conditions in a wind tunnel is extremely 
difficult. Even if it were possible to build a supersonic wind, tunnel 
which would produce an air stream large enough to take a full-size model 
and have a stagnation temperature equal to the expected flight value, 
the conditions for radiant-heat transfer would be greatly different from 
those in flight unless special care were taken to control the tempera- 
ture of the tunnel walls . 

Great simplification of experimental procedures would result if it 
were possible to interpret in terms of— corresponding flight - phenomena 
the results of tests run with considerably lower stagnation temperatures 
than those which occur in flight and with models of different size from 
the prototype. - For such tests, the models would be refrigerated or 
otherwise maintained at a temperature considerably below the normal 
stagnation temperature in order that' the test conditions would corre- 
spond to flight conditions with heat loss by radiation. The purpose of 
the present paper is to examine the similarity-conditions for flows with 
heat transfer with particular reference to the effects of model size and 
stagnation temperature and to outline a possible experimental procedure 
for obtaining the basic heat-transfer information necessary to predict 
values of the surface temperature in flight . 


SMBOLS 


F Froude number (v^/Lg) 

E heat transfer per unit area per unit“time 
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L 

M 

Nu 

Pr 

R 

Re 

T 

V 

a 

°P 

c v 

S 

k 


typical length. 


Mach number 


(y/a) 


Nusselt number (HL/k3?) 

Prandtl number ^Cp(i/k^ 
gas constant (p/pT) 

Reynolds number (pVL/p) 
absolute temperature 
typical velocity 
velocity of sound 
specific heat at constant pressure 
specific heat at constant volume 
acceleration due to gravity 

heat conductivity (gas property), that is, heat transfer per 
unit time per unit area per unit temperature gradient (fluid 
at rest) 

p pressure 

u, v, v components of velocity in rectangular coordinates 
x, y, z rectangular coordinates 


P 

7 

V- 

P 


temperature coefficient of volume expansion 

ratio of specific heat at constant pressure to specific heat 
at constant volume 

coefficient of viscosity 

density 


Subscript : 


o 


stagnation or reference conditions 
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ANALYSIS 


In examining the possibility of simulating flight conditions in a 
wind t unn el, it is of some interest to determine the effects of variations 
in the stagnation -temperature and scale of the model on the rates of heat 
transfer over the model necessary to maintain similar fields of flow, in 
particular, similar "boundary layers and similar temperature distributions 
within the boundary layers. The approach to the determination of con- 
ditions for flow similarity is similar to that of reference 2. It is 
assumed that the following conditions completely determine the field of 
flow about a body: (a) shape of the body, (b) distribution of local 

rates of heat transfer per unit area over the surface of the body, 

(c) physical properties and conditions of the fluid, (d) field of flow 
in which the body is placed, and (e) gravitational or other force field 
in which the phenomena take place. The kinds of quantities necessary 
for specifying these conditions are as follows: 

L length 

V velocity 

T absolute temperature 

H heat transfer per unit area per unit time 
k absolute heat conductivity 

pi absolute viscosity 

Cp specific heat at constant pressure 

p mass density 

p pressure 

0 temperature coefficient of volume expansion 

g acceleration due to gravity 

The number of kinds of variables required to specify the problem is 
thus 11. Since there are four fundamental types of -units - mass, length, 
time, and temperature - by the H theorem given in reference 2, seven 
Independent - nondimens ional combinations of the different- variables are 
required to specify the problem. 

The general problem of the flow of a compressible fluid with heat 
transfer can, therefore, be formulated symbolically In the following 
nondimens ional form: 
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PT, V 

pV 2 * 



l N a . I* 

\< J fflT' ^ 


pVL HL 


pc p ' Lg* ii ' kT 


= 0 


(1) 


Most of the nondlmensional parameters occurring in equation (l) are com- 
monly referred to by definite names as follows : 


C T)P 
k 

yL _ 


Lg 


= Prandtl number 
= Proude number 


■2^ = Reynolds number 
U 

— = Nusselt number 
kT 

If the fluid is assumed to obey the gas law. -2- = Constant, 

pT 



and, also, 

PT = 1 

Equation (l) may then be written 


( 


kT 

pV 2 


Pr, M, 


Re, Nu] = 0 


■) 


( 2 ) 


In accor dan ce with the definition of similar flows, all the local values 
of the parameters occurring in equation (2), as well as more directly 
descriptive ratios such as u/V, v/v, w/v, p/p 0 , ant ^- p/Po> mu3 't have 

the same values at corresponding points x/L, y/L, and z/L. If flow 
patterns are to remain simllpr as the stagnation temperature T Q and 
the size of the model L are varied, certain conditions must be fulfilled 
by the gas and the local rates of heat transfer over the body must vary 
in a specific manner. These conditions are derived from consideration 
of the various par ame ters occurring In equation (2) . 


Condition 1 - the parameter kT/piV 2 : Since T 

the local value of a 2 if y is assumed constant, 


is proportional to 
T /v 2 is proportional 
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to l/M I® and, hence, will be the same at corresponding points. The 
value of k/p must therefore be independent of* the temperature. 

Condition 2 - the parameter Pr: With the previous condition 1 

that k/| a must be constant, constancy of the Prandtl number Pr 
demands that c p be independent of the temperature . This condition is 

consistent with the previous assumption that y be constant. 

Condition 3 ” the parameter M: The value of M will be the same 

at corresponding points, provided that the free -stream velocity is varied 
as the square root of the temperature; that is, V must be proportional 

to Vv 

Condition 4 - the parameter F: Since, from condition 3, V must 

be proportional to \/t^, in order for the value of V^/Lg to remain 
unchanged, L must be proportional to T 0 (g considered constant) . 

Condition 5 - the parameter Pe : Since the velocity must vary as 

the square root of the temperature, constancy of the Reynolds number 

demands that p/p be proportional to L \/T or L pre- 
condition 6 - the parameter Nu: In order to maintain similar flow 
as L and T 0 are varied, the local rates of heat’ transfer per unit 

area from the surface H must be proportional to KT/l. 

Condition 1 and condition 2, as well as the general condition that 

= Constant, depend on intrinsic properties of the gas. The kinetic 

theory of gases composed of small hard spheres (reference 3) indicates 
that p/pT, k/p, and c p should be constant. The type of gas assumed 

in this elementary theory, therefore, satisfies the necessary conditions 
for flow similarity. In m an y respects, air behaves substantially as 
such a gas. For example, p/pT is found to be substantially constant 
over a range of temperatures from slightly above those required for 
liquefaction to those at which the effects of dissociation begin to be 
appreciable, if the pressure is not too high. On the other hand, the 
experimentally determined variation of c p , p, and k with tenperature 

is somewhat different from that predicted from the elementary theory. 

The data of reference 4 indicate a 25 "percent increase in the 
ratio k/p over the range of temperatures from -100° F to l600° F. 
Variation of k/p affects the value of kT/pV^, the- Index of the ratio 
of the flow of heat to the energy dissipated by viscous friction. 

Little 1 b known of the effect - of variations of this parameter on either 
the heat -flow pattern or the kinematic -flow pattern. It seems probable, 
however, that the primary effect of the parameter kT/pV^ is on the 
value of the surfape temperature corresponding to zero heat flow to or 
from the fluid at the boundary (equilibrium temperature). Large values 
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of k/p would correspond to low equilibrium temperatures. It has gen- 
erally been observed in heat-transfer investigations that data for widely 
different values of k /p correlate well on the basis of the Prandtl 
number. In the absence of more detailed information, it seems reasonable 
to assume that the changes in flow phenomena associated with variations 
of k/(j i will be of the same order of magnitude as the variations caused 
by corresponding changes in the Prandtl number. 

The data of reference 4 indicate that the Prandtl number for air 
varies from 0-743 at -100° F to 0.643 at l600° F. Thus, the variation 
of c.p with temperature partly compensates the variation of k/p . The 

effect on the flow pattern of a 15 -percent variation of the Prandtl number 
Is estimated to be small since the Prandtl number apparently occurs only 
to the 0.4 or 0.3 power in the relations for heat transfer given In 
reference 5- Some of this variation of the Prandtl number may be caused 
by inaccuracies in the measured values of k. The probable error in k 
is estimated in reference 4 to be 7 percent. 

In addition to the effect of a variation of Cp on the Prandtl 

number, variation of Cp Is always associated with a variation in the 

ratio of specific heats y . Variations of y cause variations in the 
pressure ratios, density ratios, and so forth, associated with given 
values of the Mach number and, therefore, can have a marked effect on 
the flow pattern. The variations of 7 become progressively more 
important above temperatures of about 700° F absolute as seen In figure 1. 
The data presented in figure 1 were derived from values Cp/R presented 

In reference 6 and the relation Cp - c T = E. 

The preceding discussion Indicates the extent to which air may be 
assumed to act as a perfect gas for the purposes of this analysis. Con- 
dition 1 and condition 2 are substantially satisfied by ai -0 provided 
that the actual air temperature in the field of flow does not vary 
sufficiently for the corresponding variations of k/p and y to have 
an appreciable effect. Variations of k/p and y will limit the range 
of values of T Q over which strict similarity is obtainable . The amount 
of variation in these ratios that Is permissible will depend in large 
measure on the types of phenomena that are of primary interest^ that Is, 
for example, whether heat transfer, boundary -layer stability, or force 
coefficients are the phenomena under investigation. 

It may be noted that the effects of variations of y are largest 
in the field of flow outside the boundary layers . The pressures about 
efficient supersonic , shapes with attached shock waves ordinarily do not 
approach the stagnation pressure very closely. The variation of y with 
tengjerature Is not expected to have a serious effect on flow similarity 
unless the temperature outside the boundary layer is somewhat higher 
than 700° F absolute . 

Conditions 3 to 6 affect mainly the free-stream and boundary con- 
ditions. As stated previously, if condition 3 and condition 4 are to be 
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satisfied simultaneously and if the acceleration due to gravity g is 
considered constant, L and T Q can no longer be varied independently 

but L must be proportional to T 0 . In most - cases of high-speed flows, 

the effects of free-convection currents of which the Froude number is 
an index are unimportant. Under these circumstances, co ndi tion 4 may 
usually be neglected and L and T 0 may be varied independently. 

Condition 6 can be interpreted as specifying the necessary local rates 
of heat transfer on the model. The variation of k with temperature, 
taken from reference 7 (the basis for the values presented in reference 4 ), 
is given in figure 2. It is seen that k varies substanti ally as a 

power of the absolute temperature ; that is, = ^ 7 -^ . Since T/T 0 

is the same at corresponding points for similar flows, the conductivity 
at all points in the field of flow is proportional to the conductivity 
at stagnation conditions . Condition 6 then reduces to the form 


E 


« 


k T 
o-'-o 

L 


0 .05 

From figure 2 it is seen that k is approximately proportional to T . 
Hence, condition 6 may be stated 


■^tunnel _ | °tunnel \ ^ ^f light 
H f light l T Of light ) L tunnel 


( 3 ) 


APPLICATION TO WIND-TUNNEL TESTING 


The preceding analysis shows the necessary conditions that must' be 
fulfilled by the free stream and at the surface of the model in order 
for the flow about the modeLto be similar to that occurring in flight. 

In most supersonic wind tunnels the absolute stagnation temperature is 
approximately equal to the ambient room, temperature; whereas in flight 
the stagnation temperature tends to be high for the higher Mach numbers. 
The actual temperature assumed by the surface of - the body in flight is 
the result of a heat balance between conduction through the boundary layer 
to the body and heat dissipation from the body by radiation. 

Simultaneous simulation of both fluid -heat-transfer and radiant- 
heat -transfer conditions in a wind tunnel is extremely difficult. In 
addition to obtaining stagnation temperatures equal to the flight values 
and cooling the tunnel walls in order to simulate flight -radiation con- 
ditions, a brief check of the similarity conditions shows that it - is 
also necessary for the linear dimensions of the model to be the same as 
those of the prototype . 
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A much, simpler procedure for simulating flight conditions in a wind 
tunnel is obtained by divorcing completely the problem of heat transmis- 
sion by radiation from that of heat transmission by surface conduction 
and then by using the similarity relations given in the preceding section 
to set up convenient tunnel operating conditions . This simplified approach 
may be illustrated by the f ollowing example . 


Consider the case of an airfoil section having sufficient internal 
heat conductivity so that the surface temperature can be considered uni- 
form. It is desired to find the characteristics of this airfoil corre- 
sponding to flight conditions at a definite value of the Mach number and 
Reynolds number. The stagnation temperature in flight is assumed to be 
known. The wind-tunnel air stream is adjusted so as to give the flight 
value of the Reynolds number and Mach number for a model of convenient 
size. The stagnation temperature of the tunnel air stream is chosen 
sufficiently low so that all heat transfer from the model by radiation 
is negligible. (If the stagnation temperature of the tunnel is equal 
to the ambient temperature in flight, the radiant -heat transfer per unit 
area in the tunnel will be approximately — i— th that in flight at a Mach 


number of k.O.) The characteristics of the airfoil and the surface 
temperature corresponding to a range of heat -transfer rates from the 
fluid to the model (various degrees of refrigeration of the model) are 
measured in the tunnel. 


The relation between the heat transfer by radiation and the surface 
temperature for flight conditions can be calculated from the Stefan- 
Boltzman law so that a curve of radiant-heat transfer from the surface 
against surface temperature can be plotted. By the use of equation (3) 
and the wind-tunnel heat -transfer data, a curve of heat transfer to the 
airfoil plotted against surface temperature is drawn for flight con- 
ditions'. The surface temperature corresponding to the intersection of 
these two curves is the surface temperature that will occur in flight; 
and the aerodynamic characteristics of the model, measured in the tunnel 
at the corresponding ratio of surface to stagnation temperatures, are 
the same as those occurring under the specified conditions in flight, 
provided that the temperatures in the field of flow are not so high that 
variations of k/p and 7 are important. 

When conditions are such that the temperatures over the surface of 
the body in flight cannot be assumed uniform, the problem of finding the 
distribution of surface temperatures in flight is much more complex. 
Nevertheless, since the laws for heat transmission by solid conduction 
and radiation are relatively simple, it seems probable that a convergent 
process of successive approximations to the surface temperature distri- 
butions can be devised. 



10 


NACA TN No. 1577 


CONCLUSIONS 


Examination of s i milarity conditions for flows of a compressible 
fluid .with heat transf er, with particular reference to variation of 
stagnation temperature and scale of the model. Indicates the following 
conclusions : 


1. In order for similar flows to be possible at various values of 
the stagnation temperature and for models of various scales, the gas 
must fulfill the following conditions: 

The ratio of the heat conductivity to the viscosity must be inde- 
pendent of the temperature, or at leas -tr must have the same value at 
corresponding points for model and flight conditions . 

The specific heat must be constant, or at least must have the same 
value at corresponding points. 

The temperature coefficient of volume expansion must be Inversely 
proportional to the absolute temperature . 

2. Maintenance of corresponding local values of the flow parameters 
Reynolds number, Mach number, Prandtl number, Froude number, and so 
forth, with variations of stagnation temperature and model scale imposes 
the following restrictions: 

The free-stream velocity must be proportional to the square root 
of the stagnation temperature. 

The kinematic viscosity must be proportional to the product of the 
model scale and the square root of the absolute stagnation temperature. 

If gravity effects are important, the scale of the model must be 
proportional, to the absolute stagnation temperature. 

3 • The range of stagnation temperatures over which similar flows 
are possible for air is limited by variations with temperature of k/n, 
the ratio of heat conductivity to viscosity, and of 7 , the ratio of 
specific heat at constant pressure to that aimonstant volume. Varia- 
tions of 7 are expected to have an Important effect on flow similarity 
only If the temperature outside the boundary layer exceeds 700° F 
absolute . 

4. Maintenance of corresponding temperature distributions over the 
surface demands that the following relation between the local rates of 
heat transfer at the surface for air be satisfied: 


H tunnel 
H f light 



^flight 

^tunnel 
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where 


H heat transfer per unit time per unit area 

T q stagnation temperature 

L typical length 

5* Subject to the lim i tations of conclusion 3, Important simplifi- 
cations of the experimental procedure for tests simulating flight at 
high Mach numbers are made possible by divorcing the fluid -heat-conduction 
problem from the heat -radiation problem and interpreting the wind-tunnel 
results by means of the foregoing similarity relations. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., December 19, 19^7 
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Figure 1.- Variation of ratio of speoifio heats Y with absolute temperature, 
extrapolated to zero pressure, calculated from the data and methods of 
reference 6. 
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